From Completely Positive Maps to the Quantum Markovian Semigroup Master 
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A central problem in the theory of the dynamics of open quantum systems is the derivation of 
a rigorous and computationally tractable master equation for the reduced system density matrix. 
Most generally, the evolution of an open quantum system is described by a completely positive 
linear map. We show how to derive a completely positive Markovian master equation (the Lindblad 
equation) from such a map by a coarse graining procedure. We provide a novel and explicit recipe for 
calculating the coefficients of the master equation, using perturbation theory in the weak-coupling 
limit. The only parameter external to our theory is the coarse-graining time-scale. We illustrate 
the method by explicitly deriving the master equation for the spin-boson model. The results are 
evaluated for the exactly solvable case of pure dephasing, and an excellent agreement is found 
within the timescale where the Markovian approximation is expected to be valid. The method can 
be extended in principle to include non-Markovian effects. 



I. INTRODUCTION 

The problem of the emergence of irreversible quantum 
dynamics from closed-system, unitary dynamics has oc- 
cupied the attention of many researchers since the birth 
of quantum mechanics ||l|-^ . It is generally believed that 
an acceptable solution is to view every quantum system 
as coupled to an environment, i.e., true quantum sys- 
tems are always "open" . The action of the environment 
is to perform measurements on the system, thus estab- 
lishing a preferred "pointer basis" and leading to deco- 
herence Within a Hamiltonian framework, recipes 

for deriving the associated reduced-dynamics were known 
since the early nineteen-sixties, starting with the Zwanzig 
projection technique: one writes down the Heisenberg 
equation of motion for the combined system-environment 
state, and then projects out the system by tracing out 
the environment degrees of freedom This yields 

an integro-differential equation involving an environment 
memory kernel, which must be subjected to approxima- 
tions in order to become useful. The two main techniques 
available are the derivation of master equations by the 
use of the Born-Markov approximation g , or a represen- 
tation in terms of path integrals and influence functionals 

. Unfortunately, in the former approach it is often un- 
clear whether or not complete positivity is preserved in 
the sequence of approximations one makes [|lO| , while in 
the latter approach one must make a semi-classical ap- 
proximation in order to obtain a tractable theory | [Tl| | . 

Complete positivity in reduced dynamics is the very 
common-sensical idea^ that the open-system dynamics 
must preserve the positivity of a system's density ma- 
trix (a necessary condition for the probability interpreta- 



tion to hold) in the presence of any other non-interacting 
system. Building on this notion, two seminal contribu- 
tions have been made. Kraus established an "operator- 
sum representation" which describes the most general 
completely-positive linear map on the density matrix of 
a quantum system 0. This is a formal representation 
of the dynamics, which has been used profit ably in the 
quantum information processing community |15[| , but is 
impractical to use for dynamics calculations. To address 
this, Lindblad has derived the most general completely- 
positive Markovian semigroup master equation for the 
dynamics of the density matrix . This master equa- 
tion can be integrated and solved to provide the time- 
development of the system density matrix. Both of these 
results were derived on the basis of axiomatic quan- 
tum mechanics. While systematically satisfying, this ap- 
proach nevertheless possesses the disadvantage that the 
resulting theories are necessarily phenomenological, in 
the sense that they contain no recipe for deriving their 
parameters from first principles Q| . 

Previous formal approaches therefore suffer from one 
of two disadvantages. Either the final equations are not 
necessarily completely positive, or they contain parame- 
ters which are not derived from first principles and must 
therefore be treated as phenomenological. In this paper 
we provide a derivation of the semigroup master equa- 
tion (SME) from the Kraus operator-sum representa- 
tion (OSR) which overcomes both of the above problems. 
Thus, the SME we derive is completely-positive (i.e., it 
is of Lindblad form) , while we can also provide a recipe 
for calculating the parameters that appear in the equa- 
tion. Our technique involves a coarse-graining procedure 
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which leaves us with just one phenomenological parame- 
ter: the coarse-graining time-scale r. The work presented 
here is a continuation and generalization of Jl7t , where 
the derivation of the SME from the OSR was provided 
for the first time. Here we verify the validity of the ap- 



proach in 1 17 1, by proving now that the resulting SME 
is, as required, completely positive (Section |l|). Further- 
more, we greatly expand the utility of the derivation by 
now also showing explicitly how to cal cula te the param- 
eters that appear in the SME (Section III). The method 
uses a perturbative expansion in the system-bath cou- 
pling strength. We apply our formalism to the simple 
example of a collection of spins coupled to a boson bath, 
and compare the result to the exact solution (Section IV). 
We conclude with an overview and assessment of possible 
extensions (Section^. 



II. FROM THE OPERATOR SUM 
REPRESENTATION OF REDUCED DYNAMICS 
TO THE SEMIGROUP MASTER EQUATION 

A. Brief Review of the Operator Sum 
Representation 

The dynamics of a quantum system S coupled to a 
bath B, which together form a closed system, evolves 
unitarily under the combined system-bath Hamiltonian 



H 



SB 



Hc«)I, 



H, 



(2.1) 



Here H5, Hs and H/ are, respectively, the system, bath 
and interaction Hamiltonians, and I is the identity op- 
erator. Assuming that S and B are initially decoupled, 
so that the total initial density matrix is a tensor prod- 
uct of the system and bath density matrices {p and pB 
respectively), the system dynamics are described by the 
reduced density matrix: 

p(0) ^p(t) -TrB[U(p®pB)Ut]. (2.2) 

Here Tr^ is the partial trace over the bath and 



U = exp(--H5i3i). (2.3) 

h 

By using a spectral decomposition for the bath, ps ^ 
Y.^ PmIm)(mI (where = 1), and introducing a ba- 

sis for the A^-dimensional system Hilbert space 

Ti., this can be rewritten in the OSR as |p3: 



K 



p(t)=^A,(t)p(0)At(t), 



(2.4) 



i=0 



where the Kraus operators {A.^} have matrix elements 
given by [||: 

[A.]„„(i) = ^^{m\{p\\J[t)\v)\n) ; i = {p,v). (2.5) 



K = Ng, where Nb is the number of bath degrees of 
freedom. Also, by unitarity of U, one derives the nor- 
malization condition 



K 



(2.6) 



which guarantees preservation of the trace of p: 
Tr[p(t)] = TrE,A,p(0)Al] = Tr[p(0) A,] - 

Tr[p(0)].The Kraus operators belong to the Hilbert- 
Schmidt space A{H) (itself a Hilbert space) of bounded 
operators acting on the system Hilbert space, and are 
represented by TV x matrices, just like p. 



B. 



Fixed-Basis Form of the Operator Sum 
Representation 



While the OSR evolution equation, Eq. (2^), is per- 
fectly general, it presents two major difficulties: (i) It is 
an evolution equation, rather than a differential equation, 
which expresses p{t) in terms of the initial condition and 
time-dependent operators. Calculating these is equiva- 
lent to diagonalizing the entire system-bath Hamiltonian. 
This is impractical in all but a very few exactly solvable 
models, (ii) It is not clear how to separate out the unitary 
evolution of the system from the possibly non-unitary 
one, which occurs from the coupling of the system to the 
bath and leads to decoherence. The reason is that in 
general, each Kraus operator will contain a contribution 
from both the unitary and the non-unitary components 
of the evolution. When one makes the assumption of 
Markovian dynamics, howe ver, a s in th e semigroup mas- 
ter equation (SME) [Eqs. (|2J^) , (^l3|) below], both of 
these problems are solved, i.e., one obtains a differential 
equation in which there is an explicit separation between 
terms leading to unitary and to non-unitary evolution. 
This provided the motivation in p7[ | to develop an al- 
ternative representation of the OSR in a form which ap- 
proaches the form of the SME as much as is possible, 
without yet making any Markovian assumption. We pro- 
vide only the main steps of this derivation here, and refer 
the interested reader to for full details. 

It is convenient for this purpose to introduce a fixed 
operator basis for A{H). Let {Ka}aLo, with Kq — I, be 
such a basis, so that the expansion of the Kraus operators 
is given by: 



M 



6„(t)KQ. 



(2.7) 



Under this expansion, the OSR evolution equation. 



Eq. (2.4), becomes 



M 



Pit) = J2 Xa/3(i)K„p(0)Ktj, 

a, (3=0 



(2.8) 
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where Xa/s (t) is the matrix with elements 



K 



Ca/3(i) = ^fc^alO^I/jW- 



The matrix x is clearly Hermitian, with positive diago- 
nal elements. With some algebraic manipulation [ p^ one 
can transform Eq. (2^) into: 



dpjt) 

at 



M 

-J[S(t),p(0)] + - J2 



a, (3=1 



([K„,p(0)Kt] + [K„p(0),Kt] 
where S(t) is the hermitian operator defined by 



(2.10) 



ih 



M 



^ [XaoWK,,-XoaWKi] . (2.11) 



Q = l 



In p7t a coarse-graining procedure wa s int roduced 
which allows to transform the exact Eq. ( 2.1C ) to the 
approximate SME. For convenience we repeat and clar- 
ify the derivation here. We consider three time-scales: (i) 
the inverse of the bath density of states frequency-cutoff 
Tc, (ii) a coarse- graining time-scale t which is essentially 
the time-scale for the bath's "memory" to disappear (the 
definition will be made more precise below), and (iii) a 
system time-scale 9 which is the typical time-scale for 
changes in the system density-matrix in the frame rotat- 
ing with the system Hamiltonian. We require that 



(2.14) 



and course-grain the evolution of the system in terms of 
'r- P] = P{jT); Xap-j = XaBin) , j an integer. Fur- 
ther, rewriting the OSR Eq. ( |2.10| ) as p(t) = A(i)p(0) 



and defining L{t) through A{t) - 
indicates time-ordering) we have 



T exp 



Eq. ( 2.10| ) is the desired result: it represents a fixed-basis 
OSR evolution equation, with a strong resemblance to 
the SME, as we now detail. 



C. From the Fixed-Basis Operator Sum 
Representation Equation to the Semigroup Master 
Equation 

1. Derivation of the Semigroup Master Equation by a 
Coarse Graining Procedure 

We recall that in the semigroup approach, under the 
assumptions of (i) Markovian dynamics, (ii) initial de- 
coupling between the system and the bath, and (iii) the 
requirement of complete positivity, the system evolves 
according to the SME ||,^: 

^ = L[p{t)] ^ -^[Us,p{t)]+Lo[pit)] (2.12) 
1 

LdW)]^- J2 aa/3([Fa,p(i)Ft]-|-[F„p(0,F]3]), (2.13) 

a, (3=1 

where aaf3 is a constant positive semi-definite matrix. 
This equation bears a clear resemblance to Eq. ( ^.10|) . 
Analyzing the differences between the S ME Eq. (|2.13| ) 
and this OSR evolution equation ( ^.ICj ) allows one to 
understand the precise manner in which the semigroup 
evolution arises from the OSR evolution under the above 
mentioned three conditions. An important difference 
between these two equations is the fact that the SME 
provides a prescription for determining p(t) at all times 
given p{t') as an initi al con dition at any other time 
t > t' > 0, whereas Eq. (2.10) determines p{t) in terms 
of p{0), i.e., at the special time t = where the system 
and the bath are in a product state. 



dpjt) 
dt 



Define = j'^r^^^^ L{s)ds, with rn = t: /J l{s)ds 



= L(t)[p(t)]. 



!li{s)ds (T 



(2.15) 



^3 

^n-1 f 



T X]j=o ^ ■ -Nfixt we make the assumption that on the 
coarse-graining time-scale r, the evolution generators 
L{t) commute in the "average" sense that [l^, L^] — 
0, Vj, k. Physically, we imagine this operation as arising 
from the "resetting" of the bath density operator over 
the time-scale r. This means that t must be larger than 
any characteristic hath time-scale, and explains the re- 
quirement Tc <C T. Under this assumption, the evolu- 
tion of the system is Markovian when t ^ t: A{t) = 

Yij=o ['''^j] ■ Further, under the discretization of the 
evolution, this product form of the evolution implies that 
Pj+i = exp [tLj] [pj]- In the limit of r ^ t we expand 
this exponential, to find that 



Pj+i - Pj 



(2.16) 



This equation is simply a discretization of Eq. ( 2.15 ) un- 
der the assumption that t <ti 9, where 9 is the time-scale 
of change for the system den sity m atrix. Notice in par- 
ticular that the RHS of Eq. ( 2.16| ) contains the average 
value of L{t) ove r the interval. Now, from the OSR evo- 
lution equation ( 2.10| ), we know the explicit form of L(t) 
over the first interval from to r. Discretizing over this 
interval we find that 



pi - po 



--\ 



M 



a, 13=1 

Lo[po], 



/3J 



[K„p(0),Kt] 

(2.17) 



where 



3 



X{s)ds. 



(2.18) 



Thus, in the sense of the course graining above we have 
arrived at an exphcit form for Lq- However, deriving an 
exphcit form for Li and for higher t erms beyond this first 
interval is impossible because Eq. ( ^.10 ) gives the evolu- 
tion in terms of yo(0). Since we have made the assumption 
that the bath "resets" over the time-scale r, we expect 
the bath to interact with the system in the same manner 
over every r-length coarse-grained interval. This is equiv- 
alent to assuming that Li = Lo,Vi (which of course is the 
most trivial way of satisfying the Markovian evolution 
condition [Li,Lj] — 0,V«,j ). Then, under the natural 
identificati on o f the K's with the F's of the SME, and 
using Eq. ( 2.16 ), one is led to the well known form of the 
semigroup equation of motion: 



dt 



M 



\ J2 (Xa/3)([K„,p(i)4] 



Q,/3=l 



[K„p(i),Kt 

(2.19) 



2. Positivity of the Coefficient Matrix 

The positive s emi-d efiniteness of the coefficient ma- 
trix in Eq. ( 2.13 ) is a sufficient condition for the 
preservation of complete positivity of the system dynam- 
ics Q. Thus, to complete the identification of Eq. ( 2.19| ) 
as a Lindblad equation, it only remains to be shown that 
(Xa/?) is positive semi-definite. To do so let us show first 
that Xaf} itself is positive semi-definite, i.e., that for any 
vector c, the matrix x satisfies c^c** > 0: 

a/3 i,a.(3 
2 



E 



>o. 



(2.20) 



where we used Eq. (2.9). Next 



I r 1 

{Xaf}) = - / Xafjdt = - (Xa/3(T) - Xa/3(0)) , (2.21) 



so that we must show that y(0 ) does not spoil the 
positivity. Now, from Eqs. (|2.3|) and (2.5) Ai(0) = 



p;:(m|u(o)|z.) 



P^l5^,ui-s = Z]aio^«"(0)^a' that 



&ja(0) = yp^5ao5i,(^i,/i) (recall Kq = Is). Thus 

Xa/3(0) = ^&jq(0)&*^(0) = ^ pA<^a0^/305i,(,i,M) 
i i=(/j,!/) 

= <5qo'5/3o- (2.22) 



But in Eq. ( p^ we are concerned with {xap) only for 
«7/9 > Ij so that finally, from Eq. ( ^.2l| ), the submatrix 
(Xct/s) with a, /3 > 1 is indeed positi ve se mi-definite. The 
important conclusion is that Eq. ( 2.1£ ) is in Lindblad 
form, i.e., it preserves complete positivity. This estab- 
lishes the validity of our result for the SME, and should 
be contrasted with projection-operator type derivations 
of the master equation ||^,|l^, which do not necessarily 
satisfy the complete positivity criterion. 



3. Separating Out the Hamiltonian 



We can write Eq. ( 2.19| ) in an alternative form which 
distinguishes between the system and bath contrib ution s 
to the unitary part of the evolution. Because Eq. ( 2.1C ) 

is linear in the Xapii) matrix, one can calculate X^pii) 
for the isolated system and hence define the new terms 
which come about from the coupling of the system to the 

bath: Xapit) = Xajl*) + Xajl*)- The terms which cor- 
respond to the isolated system will then produce a nor- 
mal -{i/K )\Rs,p{t)\ Liouville term in Eq. ( |2.19| ). Thus 
Eq. ( |2.19D can be rewritten as 

+ 2 E (Xa/3)([K„,p(04] + [Kop(i),K 

ct,/3=l 

(2.23) 

Identifying (xq/s) with O ag, a n d Kg with Fq,, this is seen 
to be equivalent to Eqs.( ^.12 )-(2.13), except for the pres- 
ence of the second term in the Liouvillian. This second 
term ^S^^^^, inducing unitary dynamics on the system, 

is referred to as the Lamb shift. It explicitly describes 
the effect the bath has on the unitary part of the system 
dynamics, and "renormalizes" the system Hamilt onian . 
It is often implicitly assumed to be present in Eq. ( 2.12 ) 

In summary, we have shown in this Section how coarse- 
graining the evolution over the bath memory time-scale 
T allows one to understand the connection between the 
OSR evolution a nd th e semigroup dynamics. The im- 
portance of Eq. ( ^.10 ) lies in the fact that it allows one 
to pinpoint the exact point at which the assumption of 
Markovian dynamics is made. Furthermore, due to the 
general likeness of its form to the SME, it provides an eas- 
ily translatable connection from the non-Markovian OSR 
to the Markovian SME. Notice also that the assumption 
of Markovian dynamics introduces an arrow of time in 
the evolution of the system, through the ordering of the 
environmental states: The system evolves through time 
in the direction of successive resettings of the bath. Addi- 
tionally, it is important to note that we have shown that 
our procedure leads to an explicitly Lindbladian (com- 
pletely posit ive) form of the SME, as written in final 
form in Eq. (|2.23|). 
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Finally, we address the question of the inclusion of non- 
Mar kovian effects. The approach presented here also of- 
fers a route to a systematic inclusion of non-Markovian 
effects, i.e., higher order dynamics which include bath 
memory terms. Such a derivation of a "post-Markovian" 
master equation is a long sought-after goal of the field 
of open quantum systems. Several attempts have been 
reported, but generally the resulting equations are not 
satisfactory because complete positivity is violated [ pO[ . 
In the context of the present approach, the formal exten- 
sion to go beyond the Markovian regime can be made by 
replacing the assumption that the evolution generators Lj 
commute to first order [see text below Eq. ( p. 15 )], by a 
higher order commutator. The derivation of this commu- 
tator and the resulting post-Markovian master equation 
is left to a future publication. 



III. EXPLICIT DERIVATION OF THE 
SEMIGROUP MASTER EQUATION 
PARAMETERS 



We can now exploit the coarse-grained first-ord er (i n 
time) perturbation expansion of the OSR, Eq. ( 2.10 ), 
made in the previous section, in order to derive the ex- 
plicit form of the parameters and operators appearing in 
the resulting SME. To do so, it turns out to be most 
convenient to work in the interaction picture (IP) de- 
fined with respect to the free system and bath Hamilto- 
nians. Let us the system-bath interaction Hamiltonian 
of Eq. (2.1) in the following perfectly general form:^ 



where {S^} and {Bq,} are the system and bath operators 
respectively, and {A^} are coupling coefhcients. In the 
IP we do not have to deal directly with the free system 
and bath Hamiltonians. However, as will be seen below, 
we do recover the Lamb shift. 



A. The Interaction Picture 

Transformation to the IP is accomplished by means of 
the unitary operator 

Ut = exp i-itUs) (E) exp (-itHs) = Us «) U^. (3.2) 



Operators in the IP will be denoted using explicit time 
dependence (and where there already was a time depen- 
dence, with an / subscript). Thus: 

H/(t) = U^H/Ut = y A„S„(<) ® B„(i) (3.3) 



where 



S„(i) = U^S„Us = 5]paMOS/3 



Bait) 



UJjB^Ui 



(3.4) 
(3.5) 



with paf3{0) = QafsiO) — <5q/3. Thc density matrix 
for the system and bath combined is denoted ptat(t) 
in the Schrodinger picture, and is transformed to the 
IP by ptot,/(i) = U^ptot(i)UT. The dynamics of 
Ptot,i{t) is governed by the unitary propagator U(t) = 
exp(— itHsB)!!^, where Hsb is the full system- 
bath Hamiltonian: ptot,i{t) = U(i)ptot,/(0)Ut(t). The 
Schrodinger and interaction pictures coincide at i = so 
that ptotj(O) = ptot(O) — p{0) (8) p_b(0). It is a standard 
exercise to show that [ElJ 



U(t) = T exp 



/ H/(r)dr 
Jo 



E 

n=l 



{-^r 



U„(t) 



where 



(3.1) U„(<) 



dtri 



dtiT{H/(ti 



(3.6) 

•HKi„)}. 
(3.7) 



The Dyson time-ordered product is defined with 
respect to any set of operators Oi{tr) as 

T{Oi(ti)---0„(i„)} = OrAtr,)---OrAtrJ, whcrC 

t^^ > tr2 > ... > ^T„ • The system density matrix in 
the IP is obtained, as usual, by tracing over the bath, 
which leads to the OSR: 

K 

piit) = Trs [Piotj(t)] = y A,(i)p(0) A|(i), (3.8) 



where the Kraus operators are now defined in the IP: 

A,(i) = yfv(M|U(t)|i/). (3.9) 



^Note that {S^} and {Ba} are not assumed to be linear 
operators, and that any interaction Hamiltonian can be de- 
composed into a sum of terms acting separately on system 
and bath. Furthermore, we allow {Sq} and {B^} to be time- 
dependent. 
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Repeating the derivation of Sections [IB, [I C we thu s ob- 
tain the very same form for the SME as in Eq. ( 2.19 ), but 
now it is a SME for the interaction representation, pi{t). 
FinaUy, the transformation back to the Schrodinger pic- 
ture is accomphshed by: 



(3.10) 



B. Perturbation Theory Expansion of the Kraus 
Operators 

Our next task is to calculate the Kraus ope rato rs . W e 
do so by using the expansion for U(t) andEqs. (3. 6), (3. 9). 
We have then: 



(3.11) 



where 



Jo JQ 
Qi,...,a ^ ^ 



t pt pt 

dtn / dtn-1... J dti X 



T 



(3.12) 



and we used [Sq. (i^), (tj)] = to separate the time- 
ordering operations. aI"-* is proportional to (Aq)", so 
that in the weak-coupling case of Aq, <C 1, we can trun- 
cate the expansion at small n. 



C. First Order Case 



First we note that from Eq. ( 3.11 ), with Kq = Is : 
bioit) — y/p^(5^i/. Now, let us calculate the expression 
for n = 1. In this case there is no need to worry about 
time ordering, and we have: 

nt 

A^J^At) = -^yp^E / dhXf,Sp{h){^,\Bl,{h)\ly) 

= Y,b^a{t)K^: (3.13) 



where the second line follows using Eqs. (3. 4), ( 3. 5) th e 
third from the fixed-basis operator expansion in Eq. (2.7), 
and we defined 



1 /■* 

rf^(t) = - / dhpai3{h)qaj{h). (3.14) 



This dimensionless quantity thus depends entirely on the 
transformation to the interaction picture, i.e., it contains 
no information on the system-bath coupling, but only 
on the internal system and bath dynamics. Next, let us 
identify — Sq (the system operators) and assume 
that our basis is trace-orthogonal: 

Tr [SlS^] = S^p/N^, (3.15) 

where Na is a normalization constant. Then 

fe,a(i) = -it^ E V(M|Ba"k>C""W a > 1- 

a' a" 

(3.16) 

Using these resuhs and Xafsit) = Y.i=^,ybia{t)b*p{t) we 
can reconstruct the x matrix: Xoo(0 — ~ 
for a > 1: 

Xao(i) = E ^^»{t)b*o{t) = E V^M^MM-aC*) 
a' a" 

where we used 

(X)b ^ Tr[pBX] = ^ P^(m|X|/x), (3.18) 

defining the bath-averaged expectation value of an arbi- 
trary operator X. Finally, for both a, P > 1 



a'a"f3'P" 



(3.19) 



where in the last line we used the completeness relation 

EJ^>H = Ib. 



Now, as shown in Eqs. (|2.21| ),( p.22D 



aa/3 = iXall) = - (Xapir) - Xaf}{Q)) 
T 



(3.20) 



unless both a = /3 = (i n wh ich case (xoo) = 0). To- 
gether with the SME Eq. ( ^.19| ) we thus have all the in- 
gredients. In particular, we can calculate the ^S^ Lamb 
shift term in Eq. ( |2.19D from Eq. ( |2.1l| ): 



a 



(3.21) 



where 



6 



= ^ A„.(B„,,)sC""M (3.22) 

a' a" 

is a correlation function which contributes to the Lamb 



shift. Note that unhke in Eq. ( 2.23 ), {Sj does not con- 
tain the system Hamiltonian, as indeed it should not in 
the IP. 

As for the decoherence term, we find 



q'q"/3'/3" 



(3.23) 



Note that both the Lamb shift parameters (j)a and the de- 
coherence parameters Uap formally depend on the coarse- 
graining time T. Since r is a "dummy" differentiation 
parameter in our theory, the dependence upon it should 
disappear in an explicit calculation. We deal with this in 
the example studied in Section 



D. Second Order Case 



Expanding the Kraus operators to second order yields: 



Aq^Aq2 / dt2 / dti X 
"'o Jo 



T [S«,(ii)S„,(t2)] (Ai|T[B„,(ti)B„,(i2)] W), 



provided that the basis of system operators is closed 
under multiplication (or more generally: is trace- 
orthonormal for all products of basis elements). The 
decoherence and Lamb shift parameters can then be cal- 
culated as in the first-order case. 



E. Putting It All Together 

In conclusion, we have derived an explicit form for the 
SME of Eq. ( |2.19D (with replacing K^). To find the 
full SME for a given problem, it is necessary to: 

1. Identify the system operators {Sq,} in the interac- 
tion Hamiltonian (and recall that these operators 
must be trace-orthogonal in our formalism) . 

2. Solve for the time-dependent system an d b at h o p- 
erators in the interaction picture [Eqs. ( |3.4|) ,(3.5)], 
and thus find T from Eq. ( |3 14| ). 

3. Calculate the expectation values of the bath oper- 
ators. The results of this step will depend on the 
initial state of the bath (e.g., thermal equilibrium, 
coherent state, squeezed state, etc.). 

4. Use the results of the previous steps to calculate 
the Lamb shift term ^S^, the decoherence matrix 
{Xaf}), and finally, to write down the SME. 

Since this SME is of the Lindblad form [fol it is guar- 



^ \p^\fi^{^i\B^^B^^\v)Tf^pl'''^^'^^{t) (3. 24)anteed to preserve positivity of the density matrix. Sys- 



/3i/32,7i72 



where 



Titr^'' (t) ^^J^ dh dhT (ti)pa,^, ih)] X 

(ti)ga.72(i2)]. (3.25) 

We need to compare this expression to the expansion of 
the Kraus operators in terms of the fixed basis, Ai{t) = 
^a=o biait)Ka [Eq. (|2.7|)]. To do so we must now extend 
the fixed basis set so that it includes product terms: 

M 



a=0 



ai=0 



ai ,02—0 



(3.26) 



tematic correcti ons m ay be derived by continuing the ex- 
pansion in Eq. ( 3.11 ) to higher orders in n. 



IV. EXAMPLE: SPIN-BOSON MODEL 

A. Pure Dephasing of Multiple Qubits 

As a concrete and simple example of the procedure 
described above, we consider the form of the SME de- 
rived for a collection of independent two-level systems 
(qubits) coupled via a phase-damping (non-dissipative) 
interaction to a boson bath. The Hamiltonians in this 
spin-boson model are [p2[: 



Comparing this expansion to Eq. (3.24) we can read off 
the second order time-dependent coefficients bi-a-^a2 as 



bi-,aia2 



J2 XpMf^\B,.'B,Mr$r'''"''' 



/3l/32,7l72 



(3.27) 



i.k 



(4.1) 

(4.2) 
(4.3) 
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where UuJq are the qubit energies, are coupling coef- 



ficients, bfc and b[ are the k bath mode annihilation 
and creation operators obeying the boson commutation 



relation 



bfc,b 



ISki, and Nfe = bj,bfe is the number 



operator. Comparing to Eq. ( |3.lD we read off the sys- 
tem operators as Sq, = cr*, and the bath operators as 
Ba = bfe. Below we deal with the required modifications 
to our treatment of the indices in order to account for 
these assignments. 

We assume that the boson bath is in thermal equi- 
librium at temperature T = 1/(^5 /5) (^s is the 
Boltzman constant). Thus the bath density ma- 



trix IS = 
/i = {ni,n2,. 



-^""^ = ^E^e-Z^^HM) (Ml, where 
...} are the numbers of quanta in 



all bath modes, i?^ is the energy of the field at a 
given occupation /i, and Z{T) = Tr[exp(— /3Hb)] — 



is the canonical partition 
function. Some useful results for the average number 
of quanta in the fc**^ bath mode and the averages of the 
creation and annihilation operators are: 



iK) 



(bfcb 



B 



Sk,i = Skj 
(bIb|>B = 



1 



: {hkhi)B = 



0. (4.4) 



We now proceed to calculate the various quantities ap- 
pearing in the SME. 



1. Form of the Interaction Representation Operators 



Formally, we need to solve Eqs. (|3.4|),(p.5|) for the time- 
dependent system and bath operators. In the present 
simple example, however it is clear that since cr* com- 
mutes with the system Hamiltonian, cr* (i) = 172(0) = (T*. 
Further, it is an elementary exercise to show that bfc(t) = 
|-,^gii^fct Therefore the interaction Hamiltonian in the IP 



is: 



H/(t) (Ale-'=*bfe + Are--'=*bt ). (4.5) 



i,k 



will use a comma to separate qubit and bath variables, 
as in a — (i^, fc), when all three indices are needed. When 
one of the indices is irrelevant it will simply be dropped. 
To separate groups of indices, such as the in T^'< , we 
will use a semicolon. For example, T-'.^'j, is short for F^''' 
with a = {iz,k), /3 = {jS^k'), 7 = where k', 

j' and ^1 are irrelevant. With these preparations we are 
now ready to calculate F^'''. By comparing al{t) = a],, 
hk{t) = bfee''^"* to Eqs. ( |!^,(m) and using the cor- 
rect index convention, we have that p"^{t) = Sn'S^z and 



Qkk" {t) 



^kk' 



-e^"^-*. Therefore: 



f; 



iz,'k (^) ~ ^ 



dtiplUti)qkk" ih) 



= <5..'4fc"<5c.e'"'=*/'sinc'"=* 



^ (4.6) 



where sine (a;) = sin (x) jx. 

Before proceeding to calculate the Lamb shift and de- 
coherence parameters, we should note that in the defini- 
tion of Sq and in Eq. (3.1), each is coupled to a 
Bq, with the same index a, whereas in the present case 
each (i.e., ti^) is coupled to both Bq, (i.e., b^.) and 
BJj. Let us briefly again suppress for clarity the ij^kl 
indices of the present example. By lin earity, the required 
modification is clearly that Eq. (3T£) should be replaced 
with 

a' a" 

+ Xl,{^,\Bl,\,.)[K'"{t)y] a>l (4.7) 



(where it is assumed that the expansion coefficients of 
Sq,, the pag {t), a re real, as in the example treated here). 
Using Eqs. ( |3.17D ,( |3.22| ) this leads to Lamb shift param- 
eters of the form: 



= 2 



E 



Re 



Aq' (Bq" 



= 0, 



(4.8) 



8. Calculation of the Lamb Shift and Decoherence 
Parameters 



In the derivation of Section III C the interaction Hamil- 



tonian was expressed as a sum over the single index a, 
in order to reduce the clutter of indices to a minimum. 
How ever , as seen from the interaction Hamiltonian of 
Eq. (4.3), in reality we need more indices. In particu- 



lar, we need to clarify the indexation of F^'i'of Eq. (3.14) 
each of the indices a, (3 and 7 may now correspond to 
a qubit position (denoted i or j), a Pauli matrix index 
(denoted ^ = y, z), or a bath mode (denoted k or I). 
Qubits variables can have both position and Pauli in- 
dices, but bath variables have only a mode index. We 



since using Eq. (L4) expectation values of creation and 
annihilation operators between numb er st ates vanish. As 
for the decoherence part, using Eq. (L4) again we find 
that (Bq//B^//)b and (B^„BJj„)b vanish, so that using 
Eqs. (|3T9|),(|3i2^), the decoherence parameters are: 



(Xa/3)=r Aq.A;, [(BQ.Bt„)sF^?"(r)F^r(r)' 

a'Q"/3'/3" 



-t-(Bt„B^.)BF^^'"(r)*F;^r(r) 



or, using the results for the spin-boson case: 



(4.9) 
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i' j' ^k'V \k"l" 



= rJj^fJ^iz AfeA^*sinc^ coth 



fc 



2 

(4.10) 



Our final result for the SME in the interaction picture 
can thus be written as 



dpi{t) 
dt 



1 



^ at^ {[al,pi{t)ai] + [alpi{t),al]) 



2^ 



where 



T \ ^ 



= — > AfeAf sinc^ coth 



(4.11) 



(4.12) 



and we reintroduced h. 

As commented above, the dependence on t should dis- 
appear after all is done. Let us see how this comes about 
within a simple continuum model. If we assume that 
|Afep only depends on Uk, we can rewrite this expression 
as an integral over lo: 

1 

h" Jo 

f3huj 



Tsmc — coth ■ 

2 2 



(4.13) 



where g{uj) is the density of states. In any reasonable 
physical model the density of states has a cut-off fre- 
quency LOc- Now, rsinc^ (ujt /2) is a function that peaks 
strongly at 0, has a width l/r and its integral is finite: 
i dwTsinc^ (wr/2) = 1, i.e., 

(5(t, cj) = — Tsinc^ {u)t/2) 

TT 

lim (5(t, lo) — 5{lo) 



But T, the coarse-graining timescale, must indeed be 
large compared to the timescale of the bath Tc — 1/lOc 
[recall Eq. (2.14)], so in this limit we can perform the 
integral and we finally get: 



al' = lim g(ijj)XHuj)\Huj)* coth 



f3hio 



(4.15) 



Thus, the dependence on r has indeed disappeared. 

We now apply this result to the case of a single two- 
level atom coupled to a harmonic bath. In the case 
of phonons and electromagnetic radiation, the interac- 
tion couples to the amphtude of the oscillators: x = 



^Jn/ (2mw)(bt + b), so that |A(a;)|^ oc l/w. At the rel- 
evant low-frequency regime we can equivalently use the 
high-temperature result: 



(btb)B = (bbt>j 



kT 



(4.16) 



For a three-dimensional crystal (or radiation field) 
giuj) oc w^. Collecting terms, we see that the limit yield- 
ing a^^ is well-defined. Decoherence depends quadrati- 
cally on the coupling, and linearly on temperature and 
on the density of low- frequency phonons. 



B. Model with Dissipation 

We now generalize our model to include dissipative 
terms. On the other hand, to keep the analysis tractable, 
we will consider the case of a single spin coupled to a bo- 
son b ath. We k eep the system and bath Hamiltonians of 
Eqs. (4.1), (|4.2| ). The new interaction Hamiltonian is: 



KM 



+a^(g) (Afe_bfe + Afe_bj^) 

(where Xk+ = A^_). Transforming to the IP we find: 
(Ta{t) = crae*"«°* and bfe(i) = bfee''^'=* where a = z,±, 
and luqz — 0, woi = T^a- As above, this translates into 



diagonal p and q [recall Eqs. (3.4),( |375| )]: 

Pap{t) = S^pe'^""' Qkk' (t) = 4fc'e*"'=*, (4.17) 
and we find for Gap' 

aa/nir) = -J E KaXlp{hlhk) BT{ujoa + ujk)T{-uJop - ujk) 
fc 

+ AfcQ'Afc/3'(bfcb^)_Br(LJoa - UJk)T{-UJoi3 + LOk). 

(4.18) 

(4.14) Here Xka' is the couphng coefficient for a]^, we already 



dropped the vanishing (b[bJ,)B and (bfebfc)^ terms, and 
T{u) ^ - f e'^^'dt = e''^^/hinc(ujT/2) . (4.19) 

T Jo 

In particular, for the diagonal terms we obtain: 

aaa(T) ^ E |Afca|^(bl.bfc)B|r(cJoa + Wfc)|^ 
fc 

+ |Afe„f (bfcbt)Blr(c^oa " c^fc)!'. (4.20) 

For Uzz this yields the same result as above. For the new 
decoherence parameters a++ and a we find: 
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a++(r) = ^ ^ |Afe+p [(b|bfe)sinc^ ((wfe - wo)t/2) The off-diagonal elements decay exponentially, with a 



+ (bfcb^)sinc2((tjfc+tJo)T/2) 
a— (r) = ((bIbfc)sinc2(K+^o)T/2) 

k 

+ (bfeb^)sinc2 (K - l^o)t/2) , ) (4.21) 
or, using the integral form: 



n Jo 



a—{T) = -2 / dwg(w)|A_(tj)|" X 
11 Jo 



{h^h)6 (w + Wo) + (bbt)^ {uj ~ ujq)) (4.22) 



With the appearance of the system's unitary time-scale 
(l/[jJo) we have to redefine our coarse-graining procedure. 
We can consider two opposite limits, where either i) the 
system energy, or ii) the interaction energy is dominant. 
These two limits correspond respectively to i) the sys- 
tem's internal unitary evolution being fast (so we are 
actually coarse-graining this out as well), and ii) the sys- 
tem's unitary evolution being slow. 



1. The fast system limit 

In this case ujqt 3> 1. The 6 functions are centered at 
ujQ and at — wq, much further from zero than their width. 
The ones at —loo thus do not contribute (off-resonance), 
so, similarly to Uzz we find: 



rate r^g). = (2a -I- ) . The diagonal elements approach 
the thermal equilibrium values Pqq'^'' and pfl'^'^, where 



^ther 

Poo 

^thcr 



(4.27) 



The exponential rate of convergence of the diagonal el- 
ements, i.e., the dissipation rate, is rJT^^^g = a__ + a++. 
Within the framework of our model, both rates depend 
linearly on temperature and quadratically on the cor- 
responding coupling strengths. The important difference 
between them is the presence of a^z in the dephasing rate. 
The parameters a+-|_ and a__ depend on the bath's den- 
sity of states at wq. Dissipation therefore can be quite 
slow in a number of important cases, for example when 
there is a gap in the phonon spectrum, or when ujo is 
actually greater than the cutoff frequency. In these cases 
only much weaker multi-phonon processes cause dissipa- 
tion. 

The parameter a^^, on the other hand, depends on 
the density of low-frequency phonons. This can be small 
only in very special circumstances (e.g., superfluidity, or 
a discrete phonon density of states as would be found in a 
quantum dot ) and its vanishing indeed usually causes 
macroscopic quantum-effects. In typical situations the 
rate of dephasing will be greater than the rate of dissi- 
pation. 

The important general conclusion is the following: If 
our coarse-graining includes the (fast) system as well, 
then the density matrix rapidly decoheres into the sys- 
tem's energy eigen-basis |2^. Then, (typically slower) it 
converges into the thermalized density matrix (which is of 
course also diagonal in the system's energy eigen-basis). 
See for a more detailed discussion of these different 
regimes. 



2n 



a++ = -2 lim g{uj)\X+{uj)\^ {h^h) 



2tt 



lim g(a;)|A_(a;)|2(bbt). 



(4.23) 
(4.24) 



The off-diagonal Lindblad parameters a^^ vanish, as they 
involve the product of two S functions that are centered 
further apart than their widths. The coefficient matrix 
is diagonal in the fast system limit. If we assume for fur- 
ther simplicity that X^{lu) — A_(cj), then the diagonal 
parameters a+-|_ and a only differ in the bath expec- 
tation values at loq. It then follows that 



(4.25) 



Let us now consider briefly the resulting interaction- 
picture master equation. Using the notation a = azz, b = 
a++, c = a for simplicity, we obtain: 



^ ^ f cpii - bpoo -(2a + ^)poi 
dt \-{2a+ ^)/3io bpoo - cpn 



(4.26) 



2. The slow system limit 

In this case ojqt <C 1. We consider onl y the zeroth ap- 
proximation, i.e., set too — 0. Using Eq. ( |4.18| ) we obtain: 

k 

+XlM^khl))\T{LOk)\\ 



or, integrating out the d function again, using Eq. (4. If) 
for the present low- frequency limit, and assuming real 
Afcct for simplicity: 



27r, 



r kT lim 



huj 



Aa(w)A^(cj) = aaaf3 



a„= I ^fcriim4^A„(c.)2 



1/2 



(4.28) 



Unlike in the fast-system case, the off-diagonal elements 
of the coefficient matrix do not vanish. Instead, in the 
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slow-system limit aap is a projection, i.e., aap is an outer 
product of the vector of components {aa\ with itself. 
This allows us to write the SME using just one Lindblad 
operator: 



This course-grained interaction operator is just a Hnear 
combination of the system operators as they appear in 
the interaction Hamiltonian, but with the dependence on 
the bath degrees of freedom aheady averaged out. Using 
G, the SME becomes: 



[Gp(t),Gt]) 



(4.30) 



Diagonahzing G and transforming p into G's eigenbasis 
then leads to uncoupled equations for the components of 
the transformed p. Therefore, in the slow system limit, 
the density matrix becomes diagonal in the eigenbasis of 
the course-grained interaction Hamiltonian (i.e., G), and 
for the rate of this decoherencc wc find: 



27r. 



'dec 



rkT lim 



5(w) 



{2\l{u,) + \\{u) + \l{Lo)) 



(4.31) 



C. Comparison of the Markovian Result to Exact 
Solution of Spin-Boson Model for Pure Dephasing 

The spin-boson model is exactly solvable in the pure 
dephasing limit, and we present the detailed solution in 
Appendix The result for an initial thermal bath is that 
the time-dependence of the off-diagonal terms is propor- 
tional to e^'"'^'*' , with 



r = ^ ^ lAfepsinc^ {LOkt/2) coth ; 



(4.32) 



This exact result holds for arbitrary times t and for both 
finite and infinite baths. 

On the other hand, recal l tha t the SME result for the 
single qubit case was [Eq. (4.12)]: 



hujk 



k ^f^B^ 



(4.33) 



This is the dephasing rate for a single qubit satisfying 
the Lindblad master equation 

^ = ^022 ([cr^p, cr^] + [0-2, pa^]) , (4.34) 

whence the off-diagonal poi oc exp(— 20^^^). 

How do these results relate to one another? We have 
in the Markovian case: 



p^r oc exp U| \>^^\'Hr,u;k) coth ^ 



(4.35) 



(4.29) whereas in the exact case: 



pTr' « exp 




S{t, LOk) coth 



2kBT 



(4.36) 



While superficially the similarity between these results is 
striking, there is nevertheless a crucial difference: The 
exact solution has recurrences, since its time-dependence 
is periodic (for a finite bath) , whereas the SME result is 
a purely exponential decay. Thus they describe very dif- 
ferent behaviors. Indeed, for small t {ujkt <C 1) the exact 
result decays as exp(— t^) (quantum Zeno effect ||26| , p7|| ), 
while the Markovian result always decays as exp(— i). 
This is of course not a surprise: the Markovian result 
cannot describe the dynamics for times shorter than the 
coarse-graining time-scale, t. 

Let us now turn to see the limit in which the two so- 
lutions do agree. To prevent recurrences in the exact 
solution we once again replace the sum over modes by an 
integral, to obtain: 

pSME^exp j du:g{u)\\{u)\H{r,u)coth^^ 

huj \ 



Poi « exp 



2n 



tj2 j '^9{^)\K^)\ 2fc^j 



The only difference is the appearance of r and t in the 
widths of the 5 functions. Now, our coarse-graining pro- 
cedure was defined such that t ^ l/ij-'c [recall the discus- 



sion surrounding Eq. (4.14)], and in this limit 5{t,uj) = 
5{ijj). For times t > t, S{t,uj) — S{uj) also holds, so we 
can summarize the condition for the exact and Marko- 
vian solutions to agree as: 



t > T > 1/Wc 



(4.37) 



To illustrate this let us consider the Debye model as a 
simple example. Then: 



g{uj) cx 



to for uj < LOc 
for Lo > LOr 



As before, let the coupling coefficient A depend on w only 
due to amplitude-coupling: |A(ix')p (x uj~^. In the high- 
temperature limit coth( ^^^j ) cx UJ^^ , so that in all we 
have 



pgf ^ cx exp -CtT / dwsinc^ (wt/2) (4.38) 



p^f cx exp ( -Ct' J dwsinc^ {ujt/2) ] , (4.39) 
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where C is the temperature-dependent coupHng- 
strength, with dimensions of frequency. Figure 1 shows 
the argument of the exponential, r(t), for the exact solu- 
tion and for the SME results, corresponding to different 
values of the course-graining time-scale, r. The curves 
corresponding to the SME solutions of course are just 
straight lines, as they all describe simple exponential de- 
cays. It is clear that the SME solutions cannot account 
for the initial transition period, but for sufficiently large 
T (in units of the bath cut-off time l/wc) the SME result 
approximates the exact solut ion very well at large times, 
in accordance with Eq. ( 4.37 ). 

Let us summarize: The Markovian approximation we 
introduced gives reliable results for times greater than 
the course-graining time-scale, which in turn must be 
greater than the bath cut-off time. It does not account 
for the initial (transitional) time evolution, and it should 
be applied in cases of an infinite bath with continuous 
spectrum. 



D. The Lamb Shift 



a general, completely positive, linear map. By coarse- 
graining this evolution over a time-scale typical of the 
bath (the inverse of the bath density-of-states frequency- 
cutoff), we showed how the Lindblad equation can be 
derived, and how its coefficients can be systematically 
calculated using perturbation theory in the system-bath 
coupling strength. This resolves an important shortcom- 
ing in the theory of open quantum systems: so far no 
practical general method was known which takes as in- 
put an interaction Hamiltonian, and then produces the 
Lindblad equation together with all its coefficients. The 
complexity of our method is determined by the difficulty 
of calculating certain time-ordered integrals, which of 
course increases with higher orders of perturbation the- 
ory. In principle, this is equivalent to the calculation 
of standard Feynman diagrams, and thus the arsenal of 
techniques known in many-body physics could be em- 
ployed here as well. To test the validity of our theory, we 
compared it here to an exactly solvable model, namely, 
the spin-boson Hamiltonian with pure phase-damping. 
For times longer than the coarse-graining time, the agree- 
ment was found to be excellent already at the level of first 
order perturbation theory. 



Finally, in the exact solution for multiple qubits there 
is also a non- vanishing Lamb shift, which arises as a con- 
sequence of the Hamiltonian not commuting with itself 
at different times . The Lamb shift does vanish for a 
single qubit in the exact solution of the pure dephasing 
spin-boson model (see Appendix ^ and [^). The Lamb 
shift also vanished in our Markovian calculation. This 
discrepancy is not only due to the fact that we considered 
a single qubit: the more fundamental reason is that we 
only carried out our Markovian calculations to first order 
in perturbation theory, where time-ordering did not play 
a role. However, when we consider the multiple-qubit 
case in second-order perturbation theory (recall Section 
HID) there is a Lamb-shift. This arises because of terms 



like a\<jlh\jDk. Physically, this is a phonon-induccd, in- 
direct, exchange-interaction between the two spins. It is 
quadratic in A, linear in temperature, and acts to pull 
the spin-energies towards an average value. 



V. CONCLUSIONS 

A central task of modern condensed phase chemistry 
and physics is the quantitative description of open quan- 
tum systems. These are systems that are coupled to an 
external uncontrollable environment (bath), a coupling 
which generally leads to decoherence. In this paper we 
provided such a quantitative description, by deriving a 
practical way to calculate the coefficients in the quan- 
tum Markovian semigroup master equation (commonly 
known as the Lindblad equation). Our starting point 
was the exact Kraus operator sum representation, which 
presents the evolution of an open quantum system as 
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APPENDIX A: ANALYTICAL SOLUTION OF 
THE SPIN BOSON MODEL FOR PURE 
DEPHASING 

We present here the analytical solution of the spin- 
boson model for pure dephasing. The derivation is based 
on[|§.^ 

Starting from the interaction picture Hamiltonian: 
Hz(t) ^Y.<® [Ale-'-'=*afe + (A^^)* e-'^'^a^] (Al) 

i,k 

we want to find the system density matrix 

piit) ^ Trs [ptotj(<)] = Trs [U(i)p(0) ® pB(0)Ut(<)] 

(A2) 

where 

ft 



U(t) = T exp 



^ Jo 



H/(r)dT 
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1. Calculation of the Evolution Operator 

Note that H/(f) does not commute with itself at differ- 
ent times, which is why we need the time-ordered prod- 
uct: 



where 



alit) 



Further: 



[Ui{t),Hi{t')] = 2i<^>ime [Xl (aQ] sina;fe {t - t') 1, 



Jo Jo 



where we used the boson commutation relations 



15ki, [a/s,ai] = 0. Note that further, 

[[Ui{t),Uj{t')],llj{t")]=0. 



and 



(A3) 



dt2 sincjk {t2 -h) = 



dti / dt2 COSWk {t2 — ti) 





This means that we can use the Baker-Hausdorf formula 

cxp{A + B) = cxp(-[A,B]/2)exp(A)exp(B) (valid if so that 



sinwkt — uJkt 



1 — cos UJkt 



cot 



(A6) 



(A7) 



(A8) 



= [[A,B],B] = 0) to calculate \J{t). To do 
so note the generalization 

exp (^A„ ) = ( Yl exp ( ) ) ( ]Jexp(A„) 



fit) 



2 „t 



\n<n' 



= ^ 2ct* o-f ® Re 



dti / dt2[H/(i2),H/(<i)] 

Wfet — sin Wfet . 



(A4) 



which is valid if every second-order commutator vanishes. 
To apply this result for our case let us formally discretize 
the integrals and denote W„ = — ^H7(nAt). Then: 



= E 



UJkt — smujkt 
{huJkf 



{huJk) 



(A9) 



V{t) = T exp 



.1 



T lim exp 
At^o 



TV 



.n=0 



lim^ H exp f-i[nn,Hn'] {AtA []exp(W„Ai) 

n<n' ^ ' n 

lim n (l - li^nM {AtA n (1 - ^uAt) 
1- J E[^"'^"'HAi) 



lim 



Note that / is an operator acting just on the system, 
and is a simple phase for the case of a single qubit. Since 
the operators commute for different modes we have as 
our final simplified result for the evolution operator: 

m = e'/(*) n^^P [< ® ("^(*)^^ - 4(*)*al) 

i,k 

(AlO) 

2. Calculation of the Density Matrix 

Now recall the definition of the coherent states. These 



(I ^\ Now recall the definition of the coherent s 

~2 V^m'Hn'] {At) I exp(^ HnAt) are eigenstates of the annihilation operator: 
n<n' J n 

i) lo '^^l/ ''^*2[Hj(t2),H7(ti)] 



exp 



exp 



a\a) = a\a). 



(All) 



They are minimum-uncertainty states in a harmonic po- 
tential, etc. As is well known, 



|a)=e-H'/2y ^ 



n) 



(A12) 



In the second line we explicitly invoked time-ordering by 
using the Campbell-Hausdorf formula to deal with the 
non-commuting problem; in the subsequent lines we used where \n) are number (Fock) states. The completeness 
the re-exponentiation trick. In the final line there is no relation for the coherent states is: 
need for explicit time-ordering left, i.e., the integrals can 
be calculated as such. We find: 



^ i,k 



^ J d^a\a){ 



a\ = l 



(A13) 



where the integration is over the entire complex plane. 
They are useful in our context since they are created by 
(A5) the displacement operator 
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D (a) = exp (aa^ — a*a) = D(— a)^ 



(A14) 



acting on the vacuum state: 

D{a)\0) = \a), (A15) 
which is clearly related to U(t). We will need the result: 

D [a) D (13) ^ exp ~ °'*^ D{a + (A16) 

which is easily derived from D (a) = exp (aa^^ — a*a), 
[a, at] — 1, and the Baker-Hausdorf formula exp(v4 + 
B) = exp(-[A,i?]/2)exp(^)exp(i?) (again, valid if 
[[A,B],A]^[[A,BIB]=0). 

Now let Rik{t) = a^(i)a^ — a^(t)*afc and consider 
exp [al ® Rik (t)] : 



exp [az (g) i?] = /g tgi 



i?2 



^ (2n)! 

n=0 ^ ' 



(2n+!)! 

n=0 ^ ' 



= Is 'S) cosh i? + cr^ (g) sinh R 
= Is ® ^[D {a) + D i-a)] 

+<Jz'S)l[D {a)-D{-a)] 



+ (A17) 



This is an important result since it shows that depend- 
ing on whether the field is coupled to the qubit |0) or |1) 
state, the field acquires a different displacement. This 
is the source of the dephasing the qubits undergo, since 
when acting on a superposition state of a qubit, the qubit 
and field become entangled: 

exp [az (g R] (a|0) + 6|l))|/3) = a|0) (g) D {a) 1/3) 
+b\l)(g)D{~-a)\P) 
^e'-"^'^"'^^^^a\0) (g,\a + P) 

The evolution operator can be written as: 



n [10). 



(A18) 

Now assume that the boson bath is in thermal equilib- 
rium: 



1 



PB 



-/9Hb 



nr 

. k 

n7^exp(-/3nc.fcNfc), 



exp (^I3^hujk (^fc + ^j^ 



where the mean boson occupation number is: 



(A19) 



1 



(A20) 



As shown in ||], p. 122-3, this can be transformed into the 
coherent-state representation, with the result: 



PB = W_PB,k 

k 



where 



PB,k = 



r(Nfe) 



(fak exp — 



(A21) 



\ak){ak\. (A22) 



Now consider the system density matrix. Let Pxi,yi — 
\x)i{y\ where x,y — {0,1}. Since we are dealing with 
qubits the system density matrix is a sum of all possible 
tensor products of single qubit pure states, i.e., of terms 
of the form P{xi,yi} = Pxi.yi (g • • • (g PxN,yN- Thus it can 
be expanded as 



{xi,yi} 



(A23) 



Recall that we set out to evaluate p(t) ~ 
TiB [V{t)p{0) (g) pB(0)Ut(t)] . For simplicity let us now 
consider the case of a single qubit. It suffices to calcu- 
late the evolution of each of the four pure states \x){y\ 
separately. Thus 

=TrB [U(i)|x>(y| ®pB(0)Ut(t)] 



Tr? 



n 



'Diak) + \l){l\^Di~ak)] 



f:){y\®YlpB,'. 



Yl [|0)(0| ® {ai) + |1)(1| ® (-a,)] 



-ifHt) 



The terms in the three products match one-to-one for 
equal indices, so we can write everything as a product 
over a single index k. Using Tr(A (g) B) = Tr^ x TtB to 
rearrange the order of the trace and the products, and 
{—a) = D (a), we have: 



Px.yit) = 6x,oSy^oe'^^*^ 



ifHt) 



(gJ^Trs [D {ak)pB^kDi-ak)] 

k 

®[]TrB [D {ak)pB^kD{ak)] 

k 

+ 5x,i5y^^e'f^'^\l){Q\e-'f'^') 

® \[ Trs [D {-ak) ps.kD {-ak)] 

k 

+ <5,,i5,,ie'/(*)|l)(l|e-'^''(*) 

(gJlTrs [D{~ak)pB,kD{ak)]. 
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Consider the Tr^ terms: for |0)(0| and |1)(1| by cy- 
cling in the trace the displacement operators cancel and 
Trs [pB,k\ = 1- Thus, as expected the diagonal terms 
do not change [apart from the Lamb shift due to f{t)]. 
As for the off-diagonal terms (evaluating the trace in any 
complete basis): 



Trs [D (±2afc) pB,k 
(fi/3k exp 

(ff3k exp 



7r(Nfc) 
1 



7r(Nfc) 
Now: 



Y,{n\D{±2ak)m{f3k\ 

n 

{pk\D{±2ak)\l3k). 



(/3|D(±2a)|/3) 
= exp [± (a/3* - a* 13)] {I3\±a + f3) 
= exp [± {af3* - a*P)] x 

exp P* (±2a + /3) - ^ + \±2a + f3\^) 

= exp {-2\af ± 2 (a/3* - a*/3)) . 



Thus: 



Trs [D {±2ak) PB,k] = exp (-2|afc|^ 



d Pk exp 



exp (-2|afc|^ 



\Pk 



± 2 (afc/3^ - a^/3fc) 



[^(Nfc)exp(-4|afcp(Nfc))] 



exp 



exp 



exp 



4|afcn(Nfe) + - 



1 



.4|Afcp^^^^coth/3?iwfc/2 
(?iwfc) 



Thus decay of the off-diagonal terms goes as e 
with 



r = 4^|A 



jl-coswfei 



(A24) 



which is equivalent to the result that appeared in 
Eq.([4.32D above. 
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FIG. 1. Comparison of exact solution of the spin-boson model for single-qubit pure dephasing to the result obtained from 
the Markovian master equation. Strai g ht lin es correspond to the Markovian solution, which intersects the exact solution (thick 
line) at t = T, as seen from Eqs. (4. 38), (4. 39). The density of states of the boson bath is represented by the Debye model here. 
Data is plotted for C = 0.05 and Uc — 1- 
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